Quantum tunneling of semifiuxons 
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We consider a system of two semifiuxons of opposite polarity in a 0-7T-0 long Josephson junction, 
which classically can be in one of two degenerate states: jj, or []. When the distance a between the 
0-7T boundaries (semifluxon's centers) is a bit larger than the crossover distance a c , the system can 
switch from one state to the other due to thermal fluctuations or quantum tunneling. We map this 
problem to the dynamics of a single particle in a double well potential and estimate parameters for 
which quantum effects emerge. We also determine the classical-to-quantum crossover temperature 
as well as the tunneling rate (energy level splitting) between the states TJ. and ||. 
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I. INTRODUCTION 

7T- Josephson junctions (7r-JJs)i are intensively investi- 
gated as they promise important advantages for Joseph- 
son junction based electronics 2 *^, and, in particular, for 
JJ based qubits4. Nowadays a variety of technologies 
allow to manufacture such junctions^iLSiS. 

One can also fabricate so-called long Josephson ex- 
junctions (0-7r-LJJs)i£, i.e., LJJs some parts of which be- 
have as O-junctions and other parts as 7r-junctions. The 
most interesting fact about such junctions is that a vortex 
of supercurrent, carrying one half of the magnetic flux 
quantum <£> ps 2.07 x 10~ 15 Wb, can be formed at the 
boundaries between and tt regions. Classically, this so- 
called semifluxoniiii^ has a degenerate ground state of ei- 
ther positive or negative polarity. The difference between 
positive and negative polarity is in the direction of the 
circulation of the supercurrent and, therefore, in the di- 
rection of the resulting magnetic field. The classical prop- 
erties of semifiuxons are under intense theoretical and 
experimental investigationo 13 i 14 i 15 i 16 i 17 i 18 i 19 i 20 i 21 i 22 i 23 i 24 . 
While the classical properties of semifiuxons (at least for 
systems with few semifiuxons) are more or less under- 
stood, their quantum behavior and their possible appli- 
cations in the quantum domain still have to be studied. 

When the energy barrier separating two degenerate 
classical states is very small, the system may sponta- 
neously switch from one state to the other due to thermal 
excitation over the barrier or due to quantum tunnel- 
ing through the barrier. Thermally induced flipping of 
a single semifluxon was already observed^. The quan- 
tum tunneling in the system of two coupled semifiuxons 
was investigated theoretically by Kato and Imada™ for 
the case of a biased junction, i.e., when the ground states 
are not degenerate and escape takes place in a certain di- 
rection. In view of possible applications and fundamental 
studies, it is interesting to see how degenerate semifluxon 
systems behave when the quantum effects start to exhibit 
themselves. 

In this paper, we study the two simplest systems: 
(a) one semifluxon with degenerate states | an d J. and 



(b) two coupled, antiferromagnetically (AFM) arranged 
semifiuxons with degenerate states || and IT- I n the 
first case, we use simple arguments to show that a single 
semifluxon is always deep in the classical limit. For a 
system of two semifiuxons, we map the problem to the 
dynamics of a single particle in a double well potential 
and estimate relevant parameters for emergence of quan- 
tum effects. We also estimate the crossover temperature 
as well as the tunneling rate (energy level splitting) be- 
tween the states |j and J,f. 



II. MODEL 

We consider a long one dimensional Josephson junction 
where the Josephson phase n(x,t) is a continuous func- 
tion of the coordinate x along the LJJ and of time t. The 
dynamics of such a system is described by a Lagrangian 
L = K — U, where 



K = E, 



-2 A*t 



dx, 



(1) 



represents the kinetic energy and 



U = E.j 



,2 Ma 



+ [1 - cos( M + 6(x))] } dx, (2) 



is a potential energy. The subscripts x and t denote the 
partial derivatives with respect to coordinate and time, 
accordingly. In the above equations the three physical 
parameters are the Josephson energy per unit of junc- 
tion length Ej, the Josephson penetration depth Aj and 
the Josephson plasma frequency uo p . The function 9{x) 
describes the position of 0- and 7r-regions along the junc- 
tion. It is zero along 0-regions and is equal to tt along 
7r-regions. 

It is straightforward to derive the equations of motion 
for the Josephson phase from the Lagrangian using the 
Euler-Lagrangc prescription. Thus, on the classical level, 
one finds that the dynamics of the Josephson phase is 
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FIG. 1: Two solutions fJ,f{x) and Hi(x) corresponding to the 
two different states | and | of a single semifluxon in a 0-7T-LJ J. 
The background color shows the corresponding Josephson po- 
tential energy density U = 1 — cos(/i + 9), black corresponds 
to the valleys IA = 0, while white corresponds to the summits 
U = 2. 



described by the time-dependent sine-Gordon equation^ 



>?jlLxx - u p 2 Htt - sin[/^ + 6{x)} = 0. 



(3) 



Here, damping and bias current are absent because dis- 
sipation and driving are not included in the initial La- 
grangian. For the present discussion these terms are not 
required as we consider undriven dissipationless systems. 



III. A SINGLE SEMIFLUXON IN A 
0-7r- JUNCTION 

Let us consider an infinite 0-7r LJJ. In this case 9{x) is 
a step function 



6(x) 



0. 

7T. 



x < 0; 
x > 0. 



(4) 



Classically, the ground state of this system is a single 
semifluxon ^SiiiiiSiiS Such a semifluxon may have posi- 
tive or negative polarity that corresponds to two classical 
degenerated states | and 

Each scmifluxon's ground state can be considered as 
a string fi(x) laying in the potential profile U(x,fi) = 
l — cos[fi + 8(x)] (Josephson energy density) . This poten- 
tial profile looks like two sets of parallel valleys shifted 
relative to each other along (i direction by 7r at x = 0, 
see Fig.^ At x — > — oo the scmifluxon's string is located 
in the valley fj, = and at x — > +oo it lays in the valley 
/j = ±7r for the states f and I, respectively. This is shown 
in Fig. n by two curves corresponding to the states | and 



Up is the 



|. It is obvious that for quantum tunneling between | and 
J, states the semi-infinite half of the string x > should 
tunnel from one valley to the other over the potential 
barrier with the height (per unit of junction length) of 
the order of 2Ej. For finite LJJ, the energy barrier scales 
proportionally to the LJJ's length L, i.e., AU ~ EjL and 
the probability of tunneling is exponentially small. For 
typical parameters (see discussion in Sec. lIVDll the ther- 
mal escape exponent AU/k B T > 100Z/Aj at T = 4K 
and it becomes even larger at smaller T. The quantum 
escape exponent AU/fku n ~ 400, where ui 
eigenfrequency of a semifluxon. 

The large barrier height results from the large length of 
the LJJ. As a modification, one can consider a LJJ of fi- 
nite, rather small length L < Xj. In this case the barrier 
height is finite and approaches zero when the junction 
length L — > 0. In this limit, one can not really speak 
about a semifluxon. The solution for the phase in such 
a short junction can be found using the image technique 
(see Fig. 6 of Ref. l2lf) . It represents a fragment of an 
infinite chain of antiferromagnetically (AFM) arranged 
scmifluxons. The flux $ present in the junction is much 
smaller than <&o/2- The possibility of quantum tunnel- 
ing in such a chain or in its fragments will be studied 
elsewhere. 

For L > A j , another option for flipping from state f to 
I is the process of emitting a fluxon ff, i.e., t=l + it- This 
process takes place already on the classical leveti&Si and 
one does not need infinite energy to put the string from 
one valley to the other. However, it still requires a rather 
high energy ~ 8Ej\j to create a fluxon. Consequently, 
this scenario will not be realized cither. 

If we consider the discontinuous Josephson phase 
4>{x, t) = fj,(x, t) + 0(a;)ii instead of the continuous phase 
H, the problem outlined above remains, but the semi- 
infinite tail should tunnel between 0(+oo) = and 
0(+oo) = 2tt. 

Thus, we have shown that a single semifluxon in a LJJ 
is always in the classical limit if L > Aj. Therefore, be- 
low we consider the more complex system of two coupled 
scmifluxons, in which the barrier separating two classical 
states can be made quite small and quantum effects may 
emerge. 



IV. TWO COUPLED SEMIFLUXONS IN O-tt-0 
JUNCTION 

Let us consider an infinite 0-7T-0 LJJ. In this case the 
function 0(x) is a step function with +n and — n discon- 
tinuities situated at x — ±a/2: 



6{x) 



0. 



a 

> 2 ; 



x < 



(5) 



where a is the length of the 7r-region between 0-7T bound- 
aries. The ground state in such a junction crucially de- 
pends on a. If the distance a is smaller than a crossover 
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FIG. 2: Two solutions fifi(x) and nit(x) corresponding to the 
two different states tl and || of two semifluxons in 0-7T-0-LJJ. 
The background color shows the corresponding Josephson po- 
tential energy U = 1 — cos(/i+#), black corresponds to the val- 
leys IA — 0, while white to the summits IA = 2. The solid line 
shows the exact numerical solution of sine-Gordon Eq. 10, 
while the dashed line (almost undistinguishable from the solid 
line) shows the approximate solution JSJ for a = 1.7 'Xj for 
B = +B (state U). 



distance a c = (n/2)Xj, the ground state of the sys- 
tem is the so-called flat phase state fj,(x) = 0, while for 
a > a c two antiferromagnetically (AFM) ordered semi- 
fluxons form the ground state>i£*22i2£i Due to symmetry 
reasons, there are two possible semifluxon states fj. and 
It that have the same energy, see Fig.|21 We would like 
to calculate the tunneling probability or the energy level 
splitting due to the coupling between these two states. 

For a > a c appreciable tunneling takes place, if the 
energy barrier between the states tl and It is small. It 
is the case when the distance a is a bit larger than the 
crossover distance a c , i.e., a = a c + da, Sa <C Xj. The 
energy of the unstable flat phase state plays a role of a 
potential barrier. When the distance a — > a c , the bar- 
rier disappears and both states turn into a flat phase 
stato 18 ! 20 ! 25 . 



A. Collective coordinate 

In this subsection, we analyze the static solutions of the 
sine-Gordon equation corresponding to the || or || state 
and derive a simple analytic approximation for them in 
the limit Sa <C Aj. This analytic solution has an am- 
plitude parameter B, which can be used as a collective 
coordinate to map our problem onto the dynamics of a 
fictitious single particle in an effective potential. 

To analyze the static solutions we assume that /i does 
not depend on t, i.e., the second term in Eq. vanishes. 
In the limit da — Sa/Xj <C 1, i.e., when the states tl 
and |t are very similar to the flat phase state fi = 0, we 



solve the stationary sine-Gordon equation in the regions 
1 (x < -a/2), 2 (-a/2 < x < a/2) and 3 (x > a/2), 
assuming that the phase )jl(x) <C 1. In the regions 1 and 
3 we find Hi,s(x) = Ai, 3 e~ |x|/Aj (Ai >3 < 1), while in 
the region 2 we have ^(x) = B cos(.t / Xj) (B <C 1). By 
matching the boundary conditions a/2) = [ii(—a/2) 
and [i2(a/2) = ^(a/2), we can express in terms of 
B, namely Ai i3 = B cos(a/2Aj) exp(a/2Aj). Thus wc 
finally arrive at 



fi(x) = B 



a 
2X~j 

x 

a7 



e a/2\j e -\x\/X J 



X > 



X < 



2' 

a 
2' 



(G) 



Note, that solution © does not satisfy the condition for 
continuity of the derivative (j, x (±a/2). The mismatch in 
derivatives is rather small (~ BSa) and is beyond the 
order of our approximation. 

In Fig. H we compare the exact numerical solution with 
/j,(x) given by formula © for a = 1.7Xj (Sa « 0.13). We 
see that they are almost indistinguishable. In principle, 
Eq. JBJ) satisfies the linearized sine-Gordon equation for 
any B <^ 1. However, as we show later, only two values 
B = ±£>o minimize the energy of the system. By chang- 
ing the parameter B, we can make a smooth transition 
between the ground states tl and ||. Thus, we can use 
B as a collective coordinate. One can think about a fic- 
titious particle with an effective mass M, which moves 
along the coordinate B in an effective potential U(B) 
derived below in Sec. IIVBI 

To describe the dynamics of this particle, we let the B 
become a dynamic variable B(t). Thus, the shape of the 
solution fi(x) is fixed, but its amplitude depends on time. 
In this way B(t) approximately describes the dynamics of 
our system and, in particular, transition between states. 

One can view the derivation of the approximate so- 
lution and the following introduction of the collective 
coordinate from different angle. For a < a c , the flat 
phase state n(x) = is a stable solution of the static 
sine-Gordon equation and all linear response eigenmodes 
of the time-dependent sine-Gordon equation (j2J) are sta- 
ble. For a > a c the lowest eigenmode becomes unsta- 
ble. Therefore, within a linear approximation this mode 
would grow exponentially. Since the sine-Gordon equa- 
tion is a nonlinear equation, the mode amplitude B would 
not become arbitrarily large, but saturates at some value 
B = Bq. Basically, this represents a new solution (JSJ) 
for a > a c . Then we assume that the amplitude B of 
this mode is large in comparison with amplitudes of all 
other modes, but still small enough so that we only have 
to take into account terms up to the order B when we 
will calculate the energy in Sec. IIV Bl All other modes 
are neglected. A similar approach was used by Kato and 
Imada^. This assumption essentially means that Eq. © 
is a good approximation for /i(x) for 8a <C 1. Conse- 
quently, we use the lowest mode amplitude B as a collec- 
tive coordinate. Such an approximation is justified be- 



cause this lowest mode, according to the Sturm-Liouvillc 
theorem, has no zeros. 

Thus, instead of treating y(x,t) as a classical (quan- 
tum) field we can restrict ourselves to single particle clas- 
sical (quantum) dynamics in a one dimensional potential. 



B. Effective mass and potential 

To determine the effective mass M associated with the 
collective coordinate B(t) we substitute our solution (JSJ) 
into the kinetic energy Eq. . After integration we ob- 
tain 



K(B) = 



EjK 



H 

EjXj 



8w 



[1 + a + cos(5) + sin(a)] B 2 (t) 
, (4 + ir)B 2 (t), 



(7) 



where a = a/Xj. When we compare this result to the 
standard expression K(B) = MB 2 /2 for the kinetic en- 
ergy, we find that the inertial mass of a particle is given 
by 



M 



EjX 



1(4 



(8) 



Similarly, we find the effective potential U(B) by sub- 
stituting the solution (jfjj) into Eq. @. To calculate the 
integral in terms of elementary functions, we assume that 
B is small and expand the integrand up to terms ~ B 4 . 

For small Sa, we obtain the following expression 



U(B,Sd) ps EjXj 



1±1 B 4 - -SaB 2 + 25a 
128 2 



(9) 



One can see that for a given distance Sa, the potential 
energy U(B) is a double well potential with walls cx B 4 
and an energy barrier in the middle oc — B 2 , see Fig. [3] 
The potential energy has two minima at B = ±B$, where 



B = 



32 



Sa. 



(10) 



A similar result was obtained earlieB 2 ^ 2 ^. Thus, the 
two classically stable solutions tJ. and J,| correspond to 
B = ±Bq and can be seen on x-y, plane in Fig. [3J The 
corresponding potential U (B) is shown on U-B plane. 

After we have determined the mass M and the effective 
potential U(B) we can calculate two parameters which 
are important for the behavior of the system: the height 
AU of the energy barrier in U(B) and the frequency u>q 
associated with small harmonic oscillations around B = 
±Bq. The height of the energy barrier is given by 



AU(5a) = U(0, So) - U(B , Sa) 



-EjXjSa 2 . (11) 




FIG. 3: Schematic view on two classically stable states jj, and 
XT and the corresponding mapping to a single particle mov- 
ing along coordinate B in a one dimensional potential U(B). 
Quantum mechanically one should speak about probability 
density \ip(B)\ 2 to find the particle at different locations (area 
plot). 



The harmonic oscillator frequency will be the same for 
small oscillations around Bq and —Bq since the potential 
U(B) is symmetric. It is therefore sufficient to consider 
small oscillations around Bq. In order to find the har- 
monic oscillator frequency loq wc expand the potential 
U(B) around Bq and find 



U(B + SB)^ EjXj 



TT + 2Sa 




7T + 2 



. (12) 

where we have neglected terms proportional to Sa 3 and 
Sa 4 . Therefore, uiq is given by 



UJ 



2EjX, 
M 



'Sa 



'Sa. 



(13) 



To check our analytical expressions derived above, we 
have compared the stationary numerical solution of the 
full sine-Gordon Eq. © to the solution © with B = Bq 
from Eq. 1|10[1 . see Fig. We have also used the station- 
ary solution of the full sine-Gordon equation to calculate 
the energy difference AU (Sa) between the flat phase state 
and the state || and compared it to the energy barrier 
AU(Sa) given in Eq. i|ll|) . Furthermore, we have calcu- 
lated the lowest eigenfrequency u>o(Sa) of the state tl and 
compared it with Eq. (|13fl . These simulations were done 
for a in the range 1.57Aj . . . 2.00A,/ (Sa = . . . 0.43). For 
all three quantities we found excellent agreement between 
analytical expressions and numerical results in the limit 
Sa — > 0. Even for a = 2Xj (Sa = 0.43) the discrepancy 
between analytics and numerics is ~ 9 % for /x(0), ~ 31 % 
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for AC/ and ~ 21 % for uiq. Thus, our analytical approxi- 
mation JHJ with the collective coordinate B describes the 
classical dynamics of our system in the limit Sa <C 1 quite 
well. This suggests that we can also successfully use this 
collective coordinate approach for a quantum mechanical 
description. 

Before we introduce the Schrodinger equation for the 
system, we briefly summarize the classical dynamics of 
the system: The two stationary solutions ±Bq corre- 
spond to the | |-state and the It-state, see Fig. |21 The 
position dependence is taken into account by Eq. 10. 
These two solutions with B = ±Bq minimize the en- 
ergy of the system. For energies smaller than the barrier 
height AC/, B is restricted to one of the potential wells 
of U(B), Eq. jnj), and will move between two turning 
points. For sufficiently small energies B(t) describes har- 
monic oscillations with a frequency uig given by Eq. (|13f) . 



D. Estimation of quantum-to-classical crossover 

For a harmonic oscillator with mass M and frequency 
luq the width of the ground state is determined by 



(17) 



Our potential U(B), Eq. ©, is not harmonic, but for suf- 
ficiently large energy barriers we may approximate each 
potential well by an harmonic oscillator, see Eq. (|12H . We 
can then use Eq. (|17(1 to estimate the spread of a wave 
function in each of the potential wells. Using our expres- 
sions for the mass M (JSJ and the frequency ujq l|13|) . we 
find 



(SB 2 ) = 



1 



V4 + 7r EjXj y/sH 



(18) 



C. Schrodinger equation 

According to the classical picture, our system has two 
stable states fl and If corresponding to B = ±B , i.e., 
a single particle in one of the wells of the potential U (B) 
0. Now, we consider the problem quantum mechan- 
ically and ask what is the probability that the parti- 
cle tunnels, e.g., between the state corresponding to the 
classical positions —Bo and Bq. In quantum mechanics 
our collective coordinate becomes an operator. In "posi- 
tion" representation (or ^-representation) the stationary 
Schrodinger equation reads 

rf 82 mm] 



2M dB 2 



1>(B) = E^{B), (14) 



and 



where M is the effective mass defined by Eq. 
U(B) is the effective potential given by Eq. ©. 

To gain more insight, we measure the energy E in units 
of EjXj. Using the expression for the mass M, Eq. JSJ), 
we can write the Schrodinger equation (|14fl in the form 



4 \E J X J 



dB 2 



u(B) 



1>(B)=eil>(B), (15) 



where the energy eigenvalue e is defined by 
E/(EjXj) and 



u(B) = 



U(B) 
EjXj 



(16) 



is the scaled effective potential. There are only two di- 
mensionless parameters in the scaled Schrodinger equa- 
tion i|15fl : the distance between semifluxons Sa and the 
dimensionless factor fnu p / (EjXj). This factor plays the 
role of a scaled H and quantum effects will vanish in the 
limit Hlo p /(EjXj) — > 0. 

In section ^ we will investigate quantum tunneling 
based on the Schrodinger equation (|14fl . Before we start 
with these calculations we would like to make some sim- 
ple estimations to see whether we can expect to observe 
quantum effects in our system. 



Quantum effects are noticeable when the wave function 
in the left well overlaps with the wave function of the 
right well. This overlap should be appreciable, but not 
too large since otherwise two states will not be dis- 
tinguishable anymore. For rough estimation we take 
{SB 2 } > O.IB 2 as criterion for quantum behavior. Thus 
quantum effects will dominate if [we use Bq from Eq. ll(J|l ] 



(SB 2 
~~R 2 



frliJr, 



1 



16V4+¥ EjXj Sa 3/2 
2 



n 



= (2e)< 



2/x A L 1 
16V4+¥ V ^ ; V Cw 2 5a 3/2 

where we took into account the definitions 



> 0.1, 



A., 



27r/i d'jc' 



2ttJc 



E, 



2tt 



(19) 



(20) 



In Eq. I|20|) . fiod' is the inductance per square of the su- 
perconducting electrodes ([Iq is a permeability of vacuum, 
d! w 2Al, Xl is the London penetration depth), j c is the 
critical current density of the LJJ, C is the capacitance 
of the LJJ per unit of area and w is the LJJ's width. 

For typical parameters Xl = 100 nm, w = 1 /im 
and C = 4.1^F/cm 2 (Hypres21 technology with j c = 
100 A/cm 2 ) we get 



(SB 2 
Bl 



3 x 10~ 4 <5a~ 3/2 > 0.1. 



(21) 



Thus, quantum effects start to play a role for Sa < 0.02. 
It is interesting to note that, according to Eq. I|19|) . the 
occurrence of quantum effects does not depend on h& 

Using definitions l|20|l in terms of physical parameters 
of LJJ, we can express inertial mass M © as 



M 



(4 + ?r)u;C 
4VMo<?c 



2n 



2.4 x 10- 4 m P A 2 7 . 



(22) 
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To estimate the crossover temperature, we compare the 
barrier height AU (fTT|) with k B T, i.e., 



AU 

k B 



EjXj 8 
k 



-6a 2 



B 7T 



(23) 



For Sa = 0.01, we obtain T* as 130 mK, which is a rea- 
sonable value for observation using modern 3 He/ 4 He di- 
lution refrigerators. This value is also typical for other 
types of qubits based on JJsi 29 i 30 i 31 i 32 We would like to 
point out that T* oc w via Ej. This is natural since it is 
rather difficult to thermally activate a "heavy" vortex. 



U(B) 




V. QUANTUM TUNNELING 



Quantum mechanics tells us that if we start with a 
wave function which is localized at one of the two min- 
ima of the potential U (B) it will tunnel through the bar- 
rier and therefore also populate the other minimum. We 
could follow this picture by using a wave function which 
is localized in one of the minima as an initial condition 
and solve the time-dependent version of the Schrodingcr 
equation (|14|) numerically Obviously, in this approach 
the answer is not a single number since details will de- 
pend on the exact form of the initial wave function. 

We use a different approach which is based on the fol- 
lowing picture: Let us assume for the moment that the 
two wells of U (B) are separated by a sufficiently high en- 
ergy barrier. Then two lowest energy eigenvalues Eq and 
Ei differ by the tunnel splitting HA a = Ei — E which 
is small compared to hcjQ. The corresponding cigenfunc- 
tions are ipo(B) and tpi(B). The ground state tpo(B) is 
symmetric whereas the first excited state ipi(B) is anti- 
symmetric. We now use the sum and the difference to 
define ip± 



1 



V2 



[MB)±MB)} 



of the two eigenfunctions as an initial condition. These 
two wave functions are well localized in one of the wells 
of U(B), see Fig.0] The time evolution ofip±(B) is given 
by 



i>±(B,t) 



1 

71 
i 

71 



MB)e-^ Eot ±MB)e~^ Elt 



e-* Eat [MB)±MB)t 



-iAnil 



(24) 



As we can see, ijj±(B,t) coherently oscillates between 
tfj+(B) and ip-(B) (apart from an overall phase factor). 
The oscillation frequency is given by Ao. 

Our calculations are still correct if the energy barrier 
becomes smaller (or even disappears). However, the wave 
functions i[>+(B) and ip-(B) are not well localized in one 
of the wells of U(B) anymore. Instead they will signifi- 
cantly extend into the other well. Nevertheless, the os- 
cillation frequency between the states ip+(B) and ip-(B) 



i>-{B) 



MB) 




FIG. 4: The two lowest energy eigenvalues in a double well 
potential and the corresponding eigenfunctions. The ground 
state V'o(-B) is symmetric whereas the first excited state ipi (B) 
is anti-symmetric. The sum ijj+(B) and the difference ip-(B) 
of these two energy eigenfunctions are localized in the right 
or in the left well, respectively. 

is still given by Ao, although the interpretation that the 
oscillations between these two states describe a tunnel- 
ing process through a barrier may become questionable 
at some point. In the rest of this section we will calculate 
the energy splitting HA . 

A. Semiclassical limit 

In the semiclassical limit the tunnel splitting HAq can 
be calculated using standard methods like WKB^ 3 . or the 
instanton techniquo 3 ^ 3 ^ to obtain analytical results. In 
the latter the problem of finding HAq is essentially re- 
duced to determine the classical path which connects the 
two maxima of the inverted potential —U(B) [the min- 
ima of U (B)] and calculate the corresponding action. For 
a quartic double- well potential of the form © these cal- 
culations can be performed analytically. We find 3 ^ 



A = 8lu 



2AU 



exp 



16AIO 



(25) 



As we might have expected, Ao decreases exponentially 
with the barrier height AU. 

Using our expressions for AU (lllfl and ujq (|13|) . wc 
arrive at 



Se 



hA 

EjXj 



RiSa* 



hujp 

Ej\ 



■ exp 



rs EjXj 3 

-K 2 — da 2 

nuj„ 



(26) 
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FIG. 5: The energy splitting Se as a function of 5a for 
huj p / (Ej\j) = 2.4x 10 -3 . The gray line is the energy splitting 
according to the semiclassical approximation Eq. I|26|l . while 
symbols show 8e obtained by numerical solution of Eq. 1151 . 



where the two numerical factors K i and Ki are given by 



Ki 



32 



^7T(7T + 2) V 7r + 4 



1/4 



8.2, 



Ko = 



32^2(7r + 4) 
3(tt + 2) 



(27) 
(28) 



Since the factor EjXj / (huj p )8a 3 > appears in the expo- 
nent, the energy splitting hA is very sensitive to this 
number. Please note, that the inverse of this factor ap- 
pears in Eq. (|T$)l . 



B. Numerical results 

Our numerical calculations are based on the scaled 
Schrodingcr equation JTHJ. In Fig.[5]we have plotted the 
energy splitting 8s as a function of 8a for a fixed value of 
hu! p /(EjXj) — 2.4 x 10~ 3 which corresponds to the pa- 
rameters used in our estimations in Scc. lIVDl According 
to Fig. our semiclassical expression (|2T>|) describes the 
energy splitting reasonably well for 8 a > 0.01. For the 
parameters used in Sec. IIVDI our numerical calculations 
predict Ao/2n m 0.61 GHz whereas the semiclassical for- 
mula HJJ) gives A /2tt w 0.88 GHz. 

The two lowest eigenfunctions of the Schrodingcr 
Eq. (|15fl are shown in Fig. [5] for three values of 8a. For 
8a = two lowest states represent the ground and the 
first excited states of the particle in a B 4 potential with 
a relatively large spacing between the energy levels, see 
Fig. EJa). For 8a = 0.02, in Fig.[|Jc) one can see clearly 
that the wave functions become strongly localized in the 



potential minima. Therefore, the quantum tunneling is 
suppressed, see also Fig. 03 Finally, for the case Sa = 0.01 
shown in Fig. Efb), we have rather strong coupling and 
appreciable energy level splitting due to the wave func- 
tions overlap. The energy level splitting 8e as function 
of 3d is plotted in Fig. 

VI. CONCLUSION 



In conclusion, we have tried to map the problem of 
quantum evolution of the Josephson phase (quantum 
field theory) to the motion of a point-like particle in 
a double well potential. For the case of two coupled 
scmifluxons arranged antiferromagnetically, we have es- 
timated that the quantum effects start dominating when 
the length a of the 7r-region exceeds a c = -|Aj by less 
then 0.02A,, see Eq. (12 H . Corresponding frequency Arj 
of the wave function oscillation between the two states 
is ~ 1 GHz, which is good value to detect experimen- 
tally. The estimated classical-to-quantum crossover tem- 
perature T* ~ 130 mK (|23|l lays in a range accessible 
for modern 3 He/ 4 He dilution refrigerators and also rep- 
resents the typical crossover temperature for flux and 
fluxon qubita 29 . 30 . 31 ! 32 . 

We would like to mention that technology is advanced 
enough to fabricate huge arrays of 0-7r-j unctions carrying 
thousands of scmifluxonsi^. Thus, in the future it will be 
interesting to extend the results obtained here to larger- 
systems, e.g., to one or two dimensional fractional vortex 
crystals. 

On the other hand, it is also interesting to consider a 
single scmifluxon squeezed into a rather short junction. 
In this case, the system is very similar to a flux qubit with 
zero loop area and flux conservation does not prevent 
flipping between the states f and |. 

Not all problems can be so easily mapped to the sin- 
gle particle dynamics. For example, the experimentally 
relevant problem of quantum escape T— >J, + 1f at the 
scmifluxon's depinning current / — > probably will 
need more elaborate approaches. 
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